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We investigate the properties of scalar and pseudo-scalar mesons at finite temperature and 
quark chemical potential in the framework of the Nambu-Jona-Lasinio (NJL) model coupled to 
the Polyakov loop (PNJL model) with the aim of taking into account features of both chiral sym- 
] metry breaking and deconfinement. 

■ The mesonic correlators are obtained by solving the Schwinger-Dyson equation in the RPA ap- 

' proximation with the Hartree (mean field) quark propagator at finite temperature and density. 

^SJ I In the phase of broken chiral symmetry a narrower width for the a meson is obtained with respect 

to the NJL case; on the other hand, the pion still behaves as a Goldstone boson. 

When chiral symmetry is restored, the pion and a spectral functions tend to merge. The Mott 
temperature for the pion is also computed. 
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I. INTRODUCTION 



0^ [ Recentl;^, increasing attention has been devoted to study the modification of particles propagatingin a hot or dense 
medium 1,'^. The possible survival of bound states in the deconfincd phase of QCD [sl. \4. la. [gTItI. laLM. [To| has opened 
interesting scenarios for the identification of the relevant degrees of freedom in the vicinity of the phase transition 
. prl . [T^ . At the same time, renewed interest has arisen for the study of the p meson spectral function in a hot 
' medium [l4l [isl . since precise experimental data have now become available for this observable [20] . 

Oh, In this paper, we focus on the description of light scalar and pseudo-scalar mesons at finite temperature and quark 
Q^' chemical potential. Besides lattice calculations [2l|, [13, [H, [IS , high temperature correlators between mesonic current 
' operators can be studied, starting from the QCD lagrangian, within different theoretical schemes, like the dimensional 
reduction [2^, [2^ or the Hard Thermal Loop approximation [13, IH, |29[ . Actually both the above approaches rely on 
a separation of momentum scales which, strictly speaking, holds only in the weak coupling regime g ^ I. Hence they 
cannot tell us anything about what happens in the vicinity of the phase transition. 

On the other hand a system close to a phase transition is characterized by large correlation lengths (infinite in the 
case of a second order phase transition) . Its behaviour is mainly driven by the symmetries of the lagrangian, rather than 
by the details of the microscopic interactions. In this critical regime of temperatures and densities our investigation 
of meson properties is then performed in the framework of an effective model of QCD, namely a modified Nambu 
Jona-Lasinio model including Polyakov loop dynamics (referred to as PNJL model) [13, [Ml, [H, [H, [s^, [H, [1^ [13] ■ 
Models of the Nambu and Jona-Lasinio (NJL) type [11] have a lo ng history and have been extensively used to 
describe the dynamics and thermodynamics of the lightest hadrons [li, HO, [M [H, [H, [3, [H, [ii, [13, [ll]. Such 
schematic models offer a simple and practical illustration of the basic mechanisms that drive the spontaneous breaking 
of chiral symmetry, a key feature of QCD in its low-temperature, low-density phase. 

In first approximation the behavior of a system ruled by QCD is governed by the symmetry properties of the 
Lagrangian, namely the (approximate) global symmetry ^[/^(iV/) x SUii{Nf), which is spontaneously broken to 
SUv{Nf) and the (exact) SUc{Nc) local color symmetry. Indeed in the NJL model the mass of a constituent quark 
is directly related to the chiral condensate, which is the order parameter of the chiral phase transition and, hence, 
is non-vanishing at zero temperature and density. Here the system lives in the phase of spontaneously broken chiral 
symmetry: the strong interaction, by polarizing the vacuum and turning it into a condensate of quark-antiquark pairs, 
transforms an initially point-like quark with its small bare mass mo into a massive quasiparticle with a finite size. 
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Despite their widespread use, NJL models suffer a major shortcoming: the reduction to global (rather than local) 
colour symmetry prevents quark confinement. 

On the other hand, in a non-abelian pure gauge theory, the Polyakov loop serves as an order parameter for the 
transition from the low temperature, Zjv^ symmetric, confined phase (the active degrees of freedom being color-singlet 
states, the glueballs), to the high temperature, deconfincd phase (the active degrees of freedom being colored gluons), 
characterized by the spontaneous breaking of the Zat^ (center of SUc{Nc)) symmetry. 

With the introduction of dynamical quarks, this symmetry breaking pattern is no longer exact: nevertheless it is 
still possible to distinguish a hadronic (confined) phase from a QGP (deconfined) one. 

In the PNJL model quarks are coupled simultaneously to the chiral condensate and to the Polyakov loop: the model 
includes features of both chiral and Zjy^ symmetry breaking. The model has proven to be successful in reproducing 
lattice data concerning QCD thermodynamics [35| . The coupling to the Polyakov loop, resulting in a suppression of 
the unwanted quark contributions to the thermodynamics below the critical temperature, plays a fundamental role 
for this purpose. 

It is therefore natural to investigate the predictions of the PNJL model for what concerns mesonic properties. Since 
the "classic" NJL model lacks confinement, the a meson for example can unphysically decay into a qq pair even 
in the vacuum: indeed this process is energetically allowed and there is no mechanism which can prevent it. As a 
consequence, the a meson shows, in the NJL model, an unphysical width corresponding to this process. One of our 
goals is to check whether the coupling of quarks to the Polyakov loop is able to cure this problem, thus preventing 
the decay of the a meson into a qq pair. Accordingly, particular emphasis will be given in our work to the a spectral 
function. 

We compute the mesonic correlation functions in ring approximation (i.e. RPA, if one neglects the antisym- 
metrisation) with quark propagator evaluated at the Hartree mean field level. The properties of mesons at finite 
temperature and chemical potential are finally extracted from these correlation functions. We restrict ourselves to the 
scalar-pseudoscalar sectors and discuss the impact of the Polyakov loop on the mesonic properties and the differences 
between NJL and PNJL models. Due to the simplicity of the model where dynamical gluonic degrees of freedom are 
absent, no true mechanism of confinement is found (we will show that for the a meson the decay channel a ^ qq is 
still open also below Tc). 

Our paper is organized as follows: in Section |TT] we briefly review the main features of the PNJL model, how quarks 
are coupled to the Polyakov loop, our parameter choice and some results obtained in Ref. ^33\ which are relevant to our 
work. In Sections IIIII and llVI we address the study of correlators of current operators carrying the quantum numbers of 
physical mesons, and the corresponding mesonic spectral functions and propagators; we obtain the relevant formulas 
both in the NJL and in the PNJL cases, and discuss the main differences between the two models. Our numerical 
results concerning the mesonic masses and spectral functions are discussed in Section |Vl Particular attention is again 
focused on the NJL/PNJL comparison. Final discussions and conclusions are presented in Section IVTl 



II. THE MODEL 



A. Nambu ~ Jona - Lasinio model 



Motivated by the symmetries of QCD, we use the NJL model (see [43], |4J, |46j, |49| for review papers) for the 
description of the coupling between quarks and the chiral condensate in the scalar-pseudoscalar sector. We will use a 
two flavor model, with a degenerate mass matrix for quarks. The associated Lagrangian reads: 

JOnjl = qiil^df, - m)q + d {qqf + {qij5Tqf (1) 

In the above q = (m, d), m — diag{mu,mci), with m„ — rud = mo (we keep the isospin symmetry); finally r° (a = 
1,2,3) are SUf{2) Pauli matrices acting in flavor space. As it is well known, this Lagrangian is invariant under a 
global - and not local - color symmetry SU{Nc = 3) and lacks the conflnement feature of QCD. It also satisfles the 
chiral SUl{2) x SUii{2) symmetry if m = while m ^ implies an explicit (but small) chiral symmetry breaking 
from SUl{2) x SUb.{2) to SUf{2) which is still exact, due to the choice m„ — rud = mo. 

The parameters entering into Eq. ^ are usually flxed to reproduce the mass and decay constant of the pion as well 
as the chiral condensate. The parameters we use are given in Table IH together with the calculated physical quantities 
chosen to fix the parameters. The Hartree quark mass (or constituent quark mass) is m = 325 MeV and the pion 
decay constant and mass are obtained within a Hartree -I- RPA calculation. 
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A [GeV] 


Gi [GeV-2] 


mo [MeV] 


K^.^u)!'/' [MeV] 


.U [MeV] 


m„ [MeV] 


0.651 


5.04 


5.5 


251 


92.3 


139.3 



Table I: Parameter set for the NJL Lagrangian given in Eq. ((TJ and the pliysical quantities chosen to fix the parameters. 



B. Pure gauge sector 

In this Section, following the arguments given in [sol . Isil ] , we discuss how the deconfinement phase transition in a 
pure SU{Nc) gauge theory can be conveniently described through the introduction of an effective potential for the 
complex Polyakov loop field, which we define in the following. 

Since we want to study the SU {Nc) phase structure, first of all an appropriate order parameter has to be defined. 
For this purpose the Polyakov line 



L{x) = V exp 



/3 

dr A4 (x, t) 







(2) 



is introduced. In the above, A4 = iA'^ is the temporal component of the Euclidean gauge field (^4,^14), in which the 
strong coupling constant gs has been absorbed, V denotes path ordering and the usual notation (3 = 1/T has been 
introduced with the Boltzmann constant set to one [ks — !)• 
When the theory is regularized on the lattice, the Polyakov loop, 

l{x) = ^TrL(x), (3) 

is a color singlet under SU{Nc), but transforms non-trivially, like a field of charge one, under Zjy . Its thermal 
expectation value is then chosen as an order parameter for the deconfinement phase transition 52l, |53l,l5^. In fact, in 
the usual physical interpretation (ssl. l56j. {l{x)) is related to the change of free energy occurring when a heavy color 
source in the fundamental representation is added to the system. One has: 

=e-'3'^-^«(^'. (4) 

In the Z^v^ symmetric phase, {l{x)) = 0, implying that an infinite amount of free energy is required to add an isolated 
heavy quark to the system: in this phase color is confined. 

Phase transitions are usually characterized by large correlation lengths, i.e. much larger than the average distance 
between the elementary degrees of freedom of the system. Effective field theories then turn out to be a useful tool to 
describe a system near a phase transition. In particular, in the usual Landau-Ginzburg approach, the order parameter 
is viewed as a field variable and for the latter an effective potential is built, respecting the symmetries of the original 
lagrangian. In the case of the SU{3) gauge theory, the Polyakov line L{x) gets replaced by its gauge covariant average 
over a finite region of space, denoted as {{L{x))) [sOj. Note that {{L{x))) in general is not a SU{Nc) matrix. The 
Polyakov loop field: 

$(f)^(a(x))) = ^Tr,((L(f))) (5) 

is then introduced. 

Following [3^, [sO, [EH, we define an effective potential for the (complex) $ field, which is conveniently chosen to 
reproduce, at the mean field level, results obtained in lattice calculations. In this approximation one simply sets the 
Polyakov loop field <I>(x) equal to its expectation value $ =const., which minimizes the potential 



where 



4 



ao 


ai 


a2 


as 


bs 


bi 


6.75 


-1.95 


2.625 


-7.44 


0.75 


7.5 



Table II: Parameters for the effective potential in the pure gauge sector (Eq. (O). 



A precision fit of the coefRcients a^, bi has been performed in Ref. (35| to reproduce some pure-gauge lattice data. 
The results are reported in Table HT] These parameters have been fixed to reproduce the lattice data for both the 
expectation value of the Polyakov loop and some thermodynamic quantities f5§\ . The parameter Tq is the critical 
temperature for the deconfinement phase transition, fixed to 270 MeV according to pure gauge lattice findings. With 
the present choice of the parameters, $ and $ are never larger than one in the pure gauge sector. The lattice data 
in Ref. [svt show that for large temperatures the Polyakov loop exceed one, a value which is reached asymptotically 
from above. This feature cannot be reproduced in the absence of radiative corrections: therefore, at the mean field 
level, it is consistent to have $ and $ always smaller than one. In any case, the range of applicability of our model 
is limited to temperatures T < 2.5 (see the discussion at the end of the next section) and for these temperatures 
there is good agreement between our results and the lattice data for 

The effective potential presents the feature of a phase transition from color confinement (T < Tq, the minimum 
of the effective potential being at <& = 0) to color deconfinement (T > Tq, the minima of the effective potential 
occurring at <i> / 0) as it can be seen from Fig. [TJ The potential possesses the Z3 symmetry and one can see that, 
above Tq, it presents three minima (Z3 symmetric), showing a spontaneous symmetry breaking. 



C. Coupling between quarks and the gauge sector: the PNJL model 

In the presence of dynamical quarks the Z3 symmetry is explicitly broken. One cannot rigorously talk of a phase 
transition, but the expectation value of the Polyakov loop still serves as an indicator for the crossover between the 
phase where color confinement occurs ($ — > 0) and the one where color is deconfined ($ — > 1). 

The PNJL model attempts to describe in a simple way the two characteristic phenomena of QCD, namely decon- 
finement and chiral symmetry breaking. 

In order to describe the coupling of quarks to the chiral condensate, we start from an NJL description of quarks 
(global SUc{5) symmetric point-like interaction), coupled in a minimal way to the Polyakov loop, via the following 
Lagrangian (|35|)^: 

UNJL^qih^D'' -^o)q + Gi[iqqf + iqij5rqf] - T) , (8) 

where the covariant derivative reads = 9'' — iv4'' and = 6qA'~' (Polyakov gauge), with — —iA^. The strong 
couphng constant gs is absorbed in the definition of A^^^x) = gsA'^i^x)^ where is the gauge field {SUc{3)) and 
Aq are the Gell-Mann matrices. We notice explicitly that at T = the Polyakov loop and the quark sector decouple. 
In order to address the finite density case, it turns out to be useful to introduce the following effective Lagrangian: 

^'pNjL = ^PNJL + mi^q , (9) 

which leads to the customary grand canonical Hamiltonian. In the above the chemical potential term accounts for 
baryon number conservation which, in the grand canonical ensemble, is not imposed exactly, but only through its 
expectation value. Let us comment here the range of applicability of the PNJL model. As already stated in Ref. [35| . 
in the PNJL model the gluon dynamics is reduced to a chiral-point coupling between quarks together with a simple 
static background field representing the Polyakov loop. This picture cannot be expected to work outside a limited 
range of temperatures. At large temperatures transverse gluons are known to be thermodynamically active degrees 
of freedom: they are not taken into account in the PNJL model. Hence based on the conclusions drawn in (60j 
according to which transverse gluons start to contribute significantly for T > 2.5 Tc, we can assume that the range of 
applicability of the model is limited roughly to T < (2 — 3)Tc. 



^ We use here the original Lagrangian of Ref. [35|, with a complex Polyakov loop effective field, which implies that at fi the expectation 
values of <E> and $ are different. A different choice can be motivated 59] but we have checked that the calculations of the present work 
are not sensitive to this feature. 
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T = 0.26 GeV< To 
"Color confinement" 
{$> = — > No breaking of Z3 




T = 1 GeV > To 
"Color deconfinement" 
(<&) / — > breaking of Z3 




Figure 1: Effective potential in the pure gauge sector (Eq. (|6])) for two characteristic temperatures, below and above the critical 
temperature To. One can see three minima appearing above Tq. 
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D. Field equations 

1. Hartree approximation 



In this Section we derive the gap equation in the Hartree approximation, whose solution provides the self-consistent 
PNJL mass of the dressed quark. 

We start from the effective lagrangian given in Eq. ([9]). The imaginary time formalism is employed. One defines the 
vertices Tm, where M = {S,P}, in the scalar {Tg = I) and pseudo-scalar (Fp = i^^r"') channel. The diagrammatic 
Hartree equation reads: 




(10) 



where the thin line denotes the free propagator in the constant (we work in the mean field) background field ^4: 

S'o(p) — — —{jf — mo + ^^{jJL — lAi))'^, the thick hue the Hartree propagator S{p) — = — (^ — m + j^{ti — 

iA4))~^, the cross (x) the vertex Tm and the dot (♦) represents 201, the coupling constant in the scalar-pseudoscalar 
channel (indeed due to parity invariance only the scalar vertex contributes). 

Besides, 1 E I = ( ) is the Hartree self-energy and m = mg -I- S. The Hartree equation then reads: 




-00 



7l — ~00 



m-mo = 2GiTTT V / — (11) 

^ /a (27r)V-m + 70(/x-iA4) 



In all the above formulas, po = iujn and aj„ — (2n+ 1)t:T is the Matsubara frequency for a fermion; the trace is taken 
over color, Dirac and flavor indices. The symbol denotes the three dimensional momentum regularisation; we use 
an ultraviolet cut-off A for both the zero and the finite temperature contributions. Our choice is motivated by our 
wish to discuss mesonic properties driven by chiral symmetry considerations, a feature not well described if one only 
regularizes the T — Q part (in particular in the vector sector the Weinberg sum rule is not well satisfied). Through 
a convenient gauge transformation of the Polyakov line, the background field in Eq. (|lip can always be put in a 
diagonal form. This allows one to straightforwardly perform the sum over the Matsubara frequencies yielding (see 
also section llIIC|) : 



m - mo = 2G^Nf ^ / T^^fl - /(^p " M + - /(Ep + M - ^A^)] • (12) 
.^-^ J A [2tt) 



By introducing the modified distribution functions'^ and , here derived for Nc = 3 (with the usual notation 
/3-1/T): 

f<s>iEp) = 1 + 3 ($ + $e-/3(Bp+M)) e-/3(Sp+M) + e-3/3(£p+M) ' '-"^^^ 



the gap equation reads: 



d^p 2ra ^ 
A (2^^ 



m - mo - 2GiiV/iV, / 7:7^ — [1 - ft{Ep) - f^{Ep)] . (15) 



^ We will explicitly derive these quantities and their role in Sec. IIII CI 
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The latter is valid for any Nc providing one uses the corresponding f^'~ . Notice that Eq. (fTTj) . after computing 
the trace on Dirac and isospin indices, can be viewed as a generalization of the corresponding zero temperature and 
density NJL gap equation 

d'^p 1 



m - mo = SiGimNcNf , ,^ ^ ^ , 

after adopting the following symbolic replacements: 

p={po^p} {iujn + fi - iA4,p) 



(2^)4 



(16) 

(17) 
(18) 



2. Grand potential at finite temperature and density in Hartree approximation 



The usual techniques [4J, [61[ can be used to obtain the PNJL grand potential from the Hartree propagator (see 
HI): 



17 = f7($,$,m;T,^) = W($,$,r)- 
-2NfT f -^-^ (Trcln 



(m — mo)^ 

i + lV(^--^)/^ 



2N,Nf 



d^p 



A (27r) 



3 



+ Trc In 



1 + Le-(^''+^)/^] } 



(19) 



In the above formula Ep = \/p'^ + is the Hartree single quasi-particle energy (which includes the constituent 
quark mass). We then define and compute them for = 3: 



zt = Trc In 



= Ttc In 



l + Lte-(^^p-A')/T' 



{l + 3 ($ + ^e-^E,-^)/T^^ ^~{E,-^)/T ^ g-3(B,-p)/T| (20) 

|l + 3(^$ + ie"(^''+''^/'^) e"(^^+^)/^ + e"^(^''+'')/^| . (21) 



E. Mean field results 



The solutions of the mean field equations are obtained by minimizing the grand potential with respect to m, $ and 
$, namely (again below Nc = 3) 



dfl 





^(-62(T)l>-53$2+fe4$$') 



-6NfT 



d^p 



A ^27r)^ I 1 + 3 ($ + $e-(Sp-A*)/T) ^-(E^-t^J/T _^ ^-3{E^~t.)/T 

g-(£;p+M)/T 



-3iEp+ii)/T 







(-fe2(T)$-63$2 +&4$*^) 



-(iSp-M)/T 



A (27r)^ I 1 + 3 ($ + $e-(^p-'^)/^) e-^^p"'')/'^ + e-3(^p-'^)/^ 

g-2(i5p+M)/T I 



(22) 



(23) 
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and 

dm ^ ^ 

which coincides with the gap equation (llip . A complete discussion of the results in mean field approximation is given 
in 35]. For the purpose of this article, we only briefly discuss the result obtained in (ssj for the net quark number 
density, defined by the equation 



(25) 



T3 T3 d^i ' 

that we display in Fig. [2P. Note that an implicit /i-dependence of fl is also contained in the effective quark mass 



m and in the expectation values $ and $. Nevertheless, due to stationary equations (|22l [231 [241) . only the explicit 
dependence arising from the statistical factors has to be differentiated. 

One can see that the NJL model (corresponding to the $ ^ 1 limit of PNJL) badly fails in reproducing the lattice 
findings, while the PNJL results provide a good approximation for them. One realizes that, at a given value of T and 
/i, the NJL model always overestimates the baryon density, even if, for large temperatures, when in PNJL $ — * 1, the 
two models merge. 

On the other hand in the PNJL model below Tc (when $, $ 0) one can see from Eqs. (pn|) and (pij) that 
contributions coming from one and two (anti-)quarks are frozen, due to their coupling with (f> and $, while three 
(anti-)quark contributions are not suppressed even below Tc. This implies that, at fixed values of T and the PNJL 
value for Ug results much lower than in the NJL case. In fact all the possible contributions to the latter turn out to 
be somehow suppressed: the one- and two-quark contributions because of $ — > 0, while the thermal excitation of 
three quark clusters has a negligible Boltzmann factor. 

One would be tempted to identify these clusters of three dressed (anti-)quarks with precursors of (anti-)baryons. 
Indeed no binding for these structures is provided by the model. In any case it is encouraging that coupling the NJL 
Lagrangian with the Polyakov loop field leads to results pointing into the right direction. 

In the following Section we explore the PNJL results in the mesonic sector, investigating whether coupling the 
(anti-)quarks with the $ field constrains the dressed qq pairs to form stable colorless structures. 



III. MESONIC CORRELATORS 



In this Section, we address the central topic of our paper, i.e. the study of correlators of current operators carrying 
the quantum numbers of physical mesons. We focus our attention on two particular cases: the pseudoscalar iso- vector 
current 



and the scalar iso-scalar current: 



Jp^ix) = q{x)ij5T''q{x) (pion) 



Js{x) = q{x)q{x) - {q{x)q{x)) (sigma). 



(26) 



(27) 



These are in fact the channels of interest to study the chiral symmetry breaking-restoration pattern. In particular 
the scalar current represents the fluctuations of the order parameter. 

In terms of the above currents, the following mesonic correlation functions and their Fourier transforms are defined: 



Cr,^(q2) ^ t I d^xe^'^-^ (O |r (^rpix)j'J{0)) I O) = CP''iq')dai, 



and 



d'^xe"'-'' ( 



T 



(js(x) 4(0) 



(28) 



(29) 



In the above equations, the expectation value is taken with respect to the vacuum state and T is the time-ordered 
product. 



^ Indeed in [35l | a different regularization procedure was employed with respect to tfie cfioice adopted in tfiis paper. Namely, no cut-off was 
used for the finite T contribution to the thermodynamical potential. This choice was made in order to better reproduce lattice results 
up to temperatures T ~ 2Tc. In any case, for lower temperature this difference in the regularization is unimportant. In particular our 
qualitative discussion of the role of the field $ in mimicking confinement is independent of these details. 
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A. Schwinger — Dyson equations at T = /i = 

Here we briefly summarize the usual N JL results for the mesonic correlators [6^ [63l . [63 . [65l | , which we are going 
to generalize in Sec. (jlll Cp by including the case in which quarks propagate in the temporal background gauge field 
related to the Polyakov loop. The Schwinger - Dyson equation for the meson correlator C^^^'^ is solved in the ring 
approximation (RPA): 

c'""{q^) = n*"^(g2) + ^n^"f (2Gi)c^-^'^-^ (30) 

M' 

where the 

n^^'^' ^ J^^TT{TMS{p + q)TM'S{q)) (31) 

are the one loop polarizations and S{p) is the Hartree quark propagator. In terms of diagrams, one defines: 




and 



C 



MM 




Hence, we need the following (one loop) polarization functions: 

n^^(g') = j^-^Tr{S{p + q)S{q)). 
Thus, for example, for the pion channel: 

d'^p w? — p^ + 9^/4 



AiN,Nfh~2iN,Nfq^l2{q^) 



the loop integrals being: 



h 



d^p 1 
A (27r)4 p2 „ ^2 

d*p 1 



A (27r)4 [(p + gf)2 - m2] [p2 _ ^2 



By defining'* 



fp{q^) = -AiN.m'hiq^) 
and owing to the fact that the Hartree equation (fT6|) implies 

m — mo 



h = 



one shows that 163 



— + /p(r)— 

zGim TO^ 

+ /pw ) — — • 



2GiTO - - 

The explicit solutions of the Schwinger-Dyson equations in ring approximation then read: 
• Scalar iso-scalar sector 



CSS^^2, 



ss 



• Pseudo-scalar iso-vector sector 



C 



n^^(g2) 



G^^(g2) = n^^(q^) + n^^((72)(2Gi)C^^(g2) 
n^^(g2) 



G 



i-2Gin^^(g2) 



pion decay constant in the chiral limit [4411 . 
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B. NJL Schwinger-Dyson equations at finite T and fj, 

In order to study the problem at finite temperature and baryon density in the imaginary time formahsm (t = —it 
with r G [0, /?]), the r-ordered product of the operators replaces the usual time-ordering and all the expectation values 
are taken over the grand-canonical ensemble. 

One can decompose all the integrands, for example in I2, as a sum of partial fractions of the form 

^ (47) 



iujn — E + fl 

The sum over Matsubara frequencies is then computed by using 



\^=fiE-^^) (48) 

P ^ lUJn - E + fl 

where the Fermi - Dirac distribution function is given by: 

The integrals Ii and I2 (Eqs. ((37)) and ([38)) ) at finite temperature and density are then expressed as [H, HI, IstI . 



<Pp i-/(Ep-/i)-/(j;p + /z) 

27r)3 2Ep 

d^p 1 f{Ep + ^,) + fiEp-^I)-fiEp+g + ^i)-fiEp+g-^i) 



r ( ■ f d^P 1 

l~f(Ep-ti)-f{Ep+g + ti) 



{2Tr)^ 2Ep2Ep^q \uj + Ep-^-q + Ep oj — Ep+q — Ep 



(51) 



(these expression are implicitly taken at u u + irj to obtain retarded correlation functions). 

Then all the zero temperature results can be continued to finite temperature and density by a redefinition of Ii 
and 12- 

At (f = 0, the integral I2 reduces to: 

r'«J - "Va (2^ Epi.^-AEl) ^''^ 



so that we obtain: 



n^''(c.,0) = -8N,NfJ^^-^^^{l-f{Ep + t,)-f{Ep-^,)) (53) 

n-(.,o) = -8^^^//|f3^|^a-/(^. + M)-/(^.-M)). (54) 



It then follows: 



j2 — 4m^ 



xe(w2 -4TO2)e(4(A2 + m2) (55) 
(and of course, the real part is given by the Cauchy principal value of the integral). Hence: 
3mn'^^(cj,0) = 2NfNaLj^^m{-il2iuj)) 

^ NcNfu/ J^^Z^NicU, /i) 9(6.2 _ 4^2)0(4(^2 ^ ^2) _ ^2) (5g) 

Stt V 

3TOn^^(w,o) = 

NcNfiu;^ - 4m^) j lo^ - Am? 

(57) 



12 



with 



N{lo, ^) = (l - / (I - - / (I + ^)) . (58) 



C. PNJL Schwinger-Dyson equations at finite T and fj. 

Here we derived explicitly the expressions for the modified Fermi-Dirac distribution functions Eqs. (|13p and (|14p . 
Again, all the summation over Matsubara frequencies can be reduced to the sum of fractions like (I47p . By defining: 

F{Ep-fi + iA^) EE IF^ -A- (59) 

P lOJn — hp + fl — lA4 



one shows that: 



= f{Ep - ^ + i{Ai)ii) + f{Ep - /i + ^(^4)22) + f{Ep + 4(^4)33) (60) 

where {A/^ju are the elements of the diagonalized ^4 matrix. 
Let us write the Fermi-Dirac distribution function according to: 

where 

z+ EE Inl^l + e-^^-^"-^)) (62) 
can be viewed as a density of partition function. We then obtain 

Tr.FiEp-, + .A,) ^ -^E^ oE-p ' = dEp ' 

_ _ 1 ain(l + £te-/^(i^.-M)) __l dzt 

" P dEp " pdEp ^^"^^ 

where 2^ = In (l + Ltg-zSCBp-^)) f^^^ corresponding density of partition hmction in PNJL (already introduced in 
Eq.dini))- Hence 

($ + 2$e"'^(^p"^') e^'^^^f + e-3/3(Bp-M) 
Ti-cF{Ep- fi + iAi) = 3^^^^-^ ^e-Z'CSp-A')) e-/3(^^p-A') + e-3/3(£p-M) ' ^^^^ 

We can do the same for the F{Ep + fi — 1^4) case. 
Hence, we can define: 

and 

/. (i^.) ^ j-Tr.FiEp + , - ^4) - , (66) 

where 2: J and are the densities (pn| and (PT|) of the partition function in PNJL. 

The only changes in going from NJL to PNJL can then be summarized in the following prescriptions: 

^ ($ + 2$e"'^(^p"^)) e-^(^f + e-3/3(Bp-p) 

/(^^P - m) =^ /* (^^p) = ^ _^ 3 _^ ^g_;3(Ep-M)) e-/3(£p-p) + e-3/3(£;p-p) (^"^^ 

/(^^P + <") =^ /* (Ep) = 1 _^ 3 + e-/3(iSp+M) + e-3,3(£;p+p) ' (^^) 

Of course in the above the corresponding PNJL quark mass m, given by the Hartree equation with these modified 
distribution functions, should be used. 




Figure 3: Fermi - Dirac distribution function f{Ep — fj.) (valid for the NJL model) and the corresponding function f^{Ep) 
(valid for the PNJL one) as functions of p, for different temperatures. $ and m are taken at their mean field values. The 
upper lines refer to T = 0.3 GeV, the lower ones to T = 0.1 GeV. 



The functions f{Ep — fj.) and f^{Ep) are displayed in Fig. [3] for two different temperatures versus p, keeping $, $ 
and TO at the mean field values. For temperatures smaller than Tc (for example T = 0.1 GeV ~ Tc/2), the effect of 
the Polyakov loop turns out to be more relevant than for larger temperatures, close to Te- 
la discussing the PNJL results for the net quark density we already stressed the role of $ and !> in suppressing one 
and two (anti-)quark clusters in the confined phase. This also emerges in Fig. |3] comparing the PNJL and NJL curves 
at T = 0.1 GeV. Clearly the two models differ substantially when $, $ 0. On the contrary, as 1 they lead 

to similar results. 

We conclude this Section by stressing once more that the recipes given in Eqs. (p7|) and ([55)1 allow one to straight- 
forwardly generalize NJL results to the PNJL case. 



IV. MESON SPECTRAL FUNCTION AND PROPAGATOR 



In the rest of the paper we present our numerical results for the masses and spectral functions of the scalar (cr) and 
pseudoscalar (tt) mesons in a hot and dense environment. 

The spectral (or strength) function F^'^^ of the correlator (7*^^^ is defined according to: 

TpMMr -A — riMM/ , ■ -\ '^^ {u! + il], q) 

F iu;,g) = ^mC + g) - ^ , 2GinMM(^ + - (69) 

For the sake of simplicity, in the following we will consider only the zero momentum case: hence we will drop the 
dependence on q. One gets: 

F^iM = 2GiSTOn^"^^(^ + zr;) 

2Gi TT (1 _ 2Gi3fie n*"^ (cj)) V (2Gi3to nA"^(cj + zr;))^ ■ 
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For uj < 2m{T,ii), SJmll = hence the decay channel into a dressed qq pair is closed and the spectral function gets 
a bound state contribution expressed by a delta peak in correspondence of the mass of the meson. Indeed: 



2Gi ^q2 



-5{uj - Tom) ■ (71) 



and the meson mass tom is the solution of the equation 

1 - 2GiKen^^*^(mAf) = . (72) 

On the other hand, for uj > 2m{T,ii), StoII ^ and the meson spectral function gets a continuum contribution. 
Thus if the solution of Eq. (|7^ occurs above such a threshold, then the spectral function will still present a peak 
characterized by a width related to the decay channel M — > qq. In such a case the meson is no longer a bound, but 
simply a resonant state. If 3m 11 stays almost constant around the position of the peak, the spectral function is well 
approximated by a Lorentzian with a width given by: 

Ta/ = 2Gi!3mn*"'^(mA/) . (73) 

On the other hand, if 3m 11 varies with lu the solution of Equation (j72p and the maximum of the spectral function 
no longer coincide, the latter being typically below the former. In the following we choose to identify the mass of the 
meson with the maximum of the spectral function. 

We notice that (7*^*^ is a correlator of current operators which is the quantity investigated in lattice calculations. 
But one can al so g et useful information concerning qq scattering processes by extracting the meson propagator from 
the T-matrix [tO]. 

In the present framework it can be shown that the propagator for a meson is 

In the quasi-particle approximation, the above simplifies to: 

Dm{u.) ^ (75) 



"M 



where m\,j verifies the pole equation (I72p and the effective meson— quark coupling constant, 



an 



MM 



9Mqq ^ 



(76) 



is the residue at the pole. 



V. NUMERICAL RESULTS 



In this Section we present, in the PNJL model, our numerical results for the properties of the a and tt mesons in a 
hot and dense environment. 

The special role played by the a spectral function, embodying the correlations among the fiuctuations of the order 
parameter (the chiral condensate), was first pointed out in [39], within the NJL model. In particular, it was shown that 
in the Wigner (ordered) phase of chiral symmetry the a spectral function (which becomes approximatively degenerate 
with the TT one, due to chiral symmetry restoration) displays a pronounced peak, moving to lower frequencies and 
getting narrower as T — *■ Tc from above. The above excitations, characterizing the regime of temperatures slightly 
exceeding Tc, were then identified as soft modes, representing a precursor phenomenon of the phase transition. 

We will show in the following that the above qualitative features of the mesonic excitations are preserved, once the 
coupling with the Polyakov loop field is introduced in the NJL model. 



A. NJL vs. PNJL: Characteristic temperatures 



Before discussing the mesonic properties, we need to identify the characteristic temperatures which separate the 
different thermodynamic phases in PNJL and NJL. In order to define a "critical" temperature one would like to 



15 



PNJL 


= 0.256 GeV 


j,a-rmn ^ q 277 GeV 


Tmou = 0.272 GeV = IMT^ 


NJL 


= 0.194 GeV 


rj^a-min ^ Q GeV 


Tmou = 0.212 GeV = 1.09r<? 



Table III: Characteristic temperatures in the NJL and PNJL models at zero chemical potential. 



refer to the order parameters (vanishing in the disordered-symmetric phase and non-vanishing in the ordered-broken 
phase). The latter, as aheady pointed out, can be identified with the Polyakov loop <i> (if ruq — > oo) and with the 
chiral condensate (qq) (if 0), for the deconfinement and chiral phase transitions, respectively. Chiral symmetry 

restoration is also signalled by TOtt = rria- (or, strictly speaking, by the merging of tt and a spectral functions). 

In the present context, chiral symmetry is explicitly broken by the presence of a finite bare quark mass; nevertheless 
the discontinuity displayed by the order parameter or by its derivatives still allows one to define a critical line in the 
(/i — T) plane separating the two phases: for small values of the baryo-chemical potential the transition is known to 
be a continuous one (cross-over) but for larger densities it becomes first order. A critical temperature, identified with 
the maximum of —dm/dT (T^), is then commonly used in the literature: this can be applied both to NJL and PNJL. 
The PNJL model also displays a similar cross-over for the effective Polyakov loop $ (which was shown to occur at a 
temperature close to in Ref. l35|]). We also notice that, in the low-density limit, two-flavour lattice QCD shows a 
cross-over, occurring at the same critical temperature both for the deconfinement and the chiral transitions. 

A different choice for a common characteristic temperature is also possible, corresponding to the minimum of m„ (T) 
(typical temperature where the pion and sigma spectral functions start to merge) ; we denote it as T^-"^^"^ . 

In order to compare the NJL and PNJL results, it is useful to follow the evolution of the observables as functions of 
the temperature, expressed both in physical units (MeV) as well as rescaled in units of a characteristic temperature. 
For the latter we choose the corresponding in NJL and PNJL. Our choice is motivated in Sec. IV Dl 

These temperatures computed in PNJL and NJL (at ^ — 0), together with the Mott temperatures for the pion (the 
temperature at which the decay of a pion into a qq pair becomes energetically favourable), are quoted in Table Hill We 
remind that in the framework of PNJL a further critical temperature related to the deconfinement phase transition 
can be defined. Its value T* = 0.250 GeV corresponds to the $ crossover location. Worth noticing is that the value 
of the latter differs by only 6 MeV from . As it is evident from Table IIIIl the characteristic temperatures that we 
get in the present approach are much larger than the one which is usually quoted for the chiral/deconfinement phase 
transition in two flavour QCD {Tc ~ 173 MeV) on the basis of lattice calculations. They are also larger than the 
ones obtained in the PNJL model in Ref. j35j. This difference is due to our choice to regularize both the zero and 
the finite temperature contributions with a three-dimensional momentum cutoff, and also to the fact that we do not 
rescale the parameter Tq to a smaller value when we introduce quarks in the system. Nevertheless, for the purpose 
of the present work, the absolute value of the critical temperature is not important: the general properties of mesons 
that we discuss here are in fact independent of the specific value of Tc- 

B. Mesonic masses and spectral functions 

In Figs, m and we plot the masses of the cr and tt mesons, together with the Hartree quark mass and the Polyakov 
loop as functions of the temperature. The first evidence emerging from these figures is that the behavior of mesons in 
PNJL looks very similar to the corresponding one in NJL [H, [la [zHj III, [z3 (^^ can be seen in Fig. [7] where NJL 
and PNJL results are directly compared). 

In Fig. m at /i = 0, the tr— mass closely follows 2m below T ~ 0.25 GeV. Then the two curves decouple: the mass 
of the dressed quarks approaches its current value, while the mass of the a meson starts increasing. 

The TT— mass is small (it is a Goldstone boson if toq = 0) and approximately constant at low temperature. Then 
it starts to increase and tends to join the ct— mass above T ~ 0.25 GeV. Both tt and cr decay into qq as soon as 
rriM > 2m. This feature can be seen clearly in the lower panel of Fig. 31 where the width of the mesons is shown, 
together with their mass. Also at low temperatures, at variance with a realistic physical situation, the production of 
free qq pairs is allowed due to the non-vanishing width of the a meson. 

The two panels of Fig. |5l show the behavior of the mesonic masses as functions of temperature, for two different 
values of the chemical potential. For /i = 0.27 GeV the system undergoes a crossover from the low-temperature, 
chirally broken phase, to the high-temperature, chirally restored one, in analogy to what happens at /i = 0. As a 
consequence, the behaviour of the mesonic masses is very similar to the one shown at vanishing chemical potential 
(see Fig. the only difference being a lower critical temperature as /x is increased. 

The pattern changes, instead, at fi = 0.34 GeV, where a discontinuity in the masses (reflecting an analogous 
behavior of the chiral condensate (qq)) appears. This can be understood by observing that between /i = 0.27 and 
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Figure 4: Top: Masses of the a and tt as functions of the temperature, together with the Hartree quark mass and the Polyakov 
loop, in the PNJL model at p = 0. The threshold 2m is also plotted to show that the a mass is close to this value below 0.25 
GeV. Bottom: Same as before, but without the Polyakov loop and adding instead the width of the mesons, represented by 
error bars. 
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Figure 5: Masses of the a and tt as functions of the temperature, together with the Hartree quark mass, in the PNJL model at 
fi = 0.27 GeV (top) and fi = 0.34 GeV (bottom). The lower figure clearly displays a first order phase transition related to the 
discontinuity of the chiral condensate occurring at ~ 0.06 GeV. 
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0.34 GeV there exists a critical point [3g|, separating a crossover from a first order phase transition. 

In concluding this paragraph, we show the pion and a spectral functions in Fig. [B] Notice their progressive 
broadening as the temperature increases. Besides, they tend to merge for T > T^-"'-"'- ^ in the chirally symmetric 
phase, as expected. 

C. NJL vs. PNJL: Mesonic masses at /x = 

In Fig. [7] we show a direct comparison between NJL and PNJL results for the mesonic masses at = 0. According 
to the features discussed in the previous paragraph, the key quantity which governs the temperature evolution of the 
mesonic masses is the dressed quark mass. As it is evident from the figure, the main quantitative difference between 
the results of the two models is the shift of the critical temperature for the phase transition, which turns out to be 
higher in the PNJL model, with respect to the "classic" NJL one. From a qualitative point of view, there is a good 
agreement between the results of the two models: in both cases in fact, the a meson mass closely follows the behaviour 
of 2m for small temperatures, decreasing when one approaches the phase transition region. Above Tc, instead, 
increases with the temperature, merging with the pion mass. This behavior reflects chiral symmetry restoration, 
a feature which is correctly described by both models, and therefore not spoiled by the coupling of quarks to the 
Polyakov loop. 

D. NJL vs. PNJL: the cr spectral function 

In this paragraph we discuss the width of the cr spectral function for the process a — > qq, as a function of the reduced 
temperature y/y^-™*". As already anticipated, we find it convenient to rescale T by J|f-""" since, interestingly 
enough, the a spectral function computed in the PNJL model almost coincides with the one evaluated in NJL at 
T ~ j^^-mm ^ This can be clearly seen in Fig. [HI where the a spectral function at g = is plotted vs. frequency. 
Notice the broadening of the a spectral function in PNJL as compared to the NJL one when T > j|f-™"'^ pointing 
to a stronger production of "free" quarks in this regime. 

To quantitatively illustrate these features at zero chemical potential, in the upper panel of Fig. [S] we show the 
absolute values of the a meson width in the NJL and in the PNJL models, while in the lower panel we show the 
ratio between the widths evaluated in the two models. The width of the a evaluated in PNJL is smaller than the 
NJL one below and it is larger above jy-™™. This is an indication that in PNJL the decay channel a ^ qq 

is reduced at low temperatures with respect to the NJL case. On the contrary, above T|f-""" one can interpret the 
larger width of the quarks bound into the meson as a more efhcient deconfinement effect. In spite of the smallness of 
both absolute widths, the PNJL model entails a reduction up to ~ 40 % for Tc{PNJL), a step toward confinement: 
this can be seen in the lower panel of Fig. [51 where the relative width T a{P N J L) /T ^{N J L) is displayed as a function 
of the reduced temperature. 

We also notice that, in both channels and at /i = 0, Eq. fT^ no longer has a real solution for T — 0.396 GeV= 2.04T^ 
in NJL and T — 0.422 GeV= l.&bT^ in PNJL. These can be interpreted as the dissociation temperatures of the model ; 
the faster (in relative units) occurrence of dissociation in PNJL is in agreement with the larger width of the a meson 
at high temperatures. 

In Fig[5| (bottom) we also report the ratio of widths at a finite chemical potential, ji = 0.27 GeV. Here the 
overall situation is qualitatively similar to what happens at /i = 0. However the curve shows some peculiar features, 
overshooting one at small T/Tc (where the absolute widths are both, in any case, very small). This behavior does not 
lend itself to an immediate physical interpretation since the critical point, in the two models, differs not only for the 
value of the temperature but also of the chemical potential. Hence the "absolute" value /i = 0.27 GeV corresponds to 
different physical situations: a rescaling for the chemical potential should be performed, but it goes beyond the scope 
of the present work, since it involves a precise discussion of the phase diagram of the PNJL model compared to the 
NJL one. 



VI. CONCLUSIONS 

In the present work, we have investigated the properties of scalar and pseudo-scalar mesons at finite temperature 
and quark chemical potential in the framework of the Polyakov loop extended Nambu-Jona-Lasinio model. This 
model has proven to be particularly successful in reproducing two flavour QCD thermodynamics as obtained in lattice 
calculations ^351] : the coupling of quarks to the Polyakov loop produces a statistical suppression of the one- and two- 
quark contributions to the thermodynamics, thus remarkably improving the NJL model results at low temperatures. 

The present work was meant as a test of the PNJL model in the mesonic sector. On the one hand, it was important 
to check whether the role of pions as Goldstone bosons as well as the pion-cr degeneracy in the chirally restored phase 
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Figure 7: a and n masses as functions of the temperature, together with the Hartree quark mass in NJL (thin lines) and PNJL 
(thick lines) model (^ = 0). 



are still satisfied after coupling quarks to the Polyakov loop. On the other hand, it was interesting to investigate 
whether the coupling to the Polyakov loop can cure some problems of the "classic" NJL model description of mesons, 
such as the unphysical width of the a meson for the process a —^ qq in the chirally broken phase. 

Finally we also intended to generalize the NJL formalism in order to embody the Polyakov loop coupled to quarks. 
This turned out to be particularly useful in the mesonic sector. Indeed we have shown the important results that 
PNJL calculations can be directly deduced from NJL ones (not only for one loop calculations, but to all orders) simply 
by a redefinition of the usual Fermi - Dirac distribution function. 

Our work shows a perfect agreement between the NJL and PNJL results concerning the mesonic masses: in the high 
temperature phase, pions and a tend to merge, thus displaying the correct pattern for chiral symmetry restoration. 
In particular, the pions still survive as bound states up to Tmou — LOGT^?, and their Goldstone boson nature is still 
preserved in the chirally broken phase. 

As far as the a meson is concerned, no true confinement is observed in the model, since the unphysical width due 
to the decay into a qq pair is still present in the PNJL model. This does not come as a surprise, since no dynamical, 
self-coupled gluons are embodied in the model Lagrangian. In any case our results in PNJL on the decay width 
improve slightly the NJL ones. 
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